Intrinsic Spin Hall Effect in Non-Cubic Crystals 
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We study the dependence of the intrinsic spin Hall effect on the crystal symmetry and geometry 
of experiment. The spin current is obtained and the Hall voltage caused by the polarization of the 
electron spins is computed. The unique dependence of the effect on the crystal symmetry permits 
the choice of geometry in which the spin Hall effect can be unambiguously distinguished from the 
effects due to the orbital motion of charge carriers and due to the magnetic field generated by the 
transport current. 

PACS numbers: 72.25.-b, 72.10.-d, 71.70.Ej, 85.75.-d 



The spin Hall effect consists of opposite "Hall" cur- 
rents for charge carriers with opposite spin polarizations. 
Remarkably, no magnetic field is needed for the electrons 
with spin up and spin down to veer in opposite direc- 
tions. The effect was initially observed in GaAs strips by 
Kato et al. via Kerr signal and by Wunderlich et al. 
via polarization of light emitted at the recombination of 
spin-polarized electron- hole pairs 0]. When numbers of 
spin-up and spin-down electrons are equal, no net elec- 
tric current across the strip is generated by the spin Hall 
effect. If, however, the electrons are polarized due to, 
e.g., injection from a ferromagnet Q, then the numbers 
of electrons veering left and right are different and a Hall 
voltage is produced. Valenzuela and Tinkham used this 
effect to measure the spin Hall conductivity in the alu- 
minum strip 0]. 

Dyakonov and Perel [j| were apparently the first to 
suggest that scattering of charge carriers by unpolarized 
impurities in semiconductors can lead to the spin polar- 
ization of the sample boundaries. Their suggestion was 
based upon a similar effect in atomic physics known as 
Mott scattering [6(. When a relativistic electron passes 
at a speed v through an atom, the partially unscreened 
electric field of the nucleus, E, creates the magnetic field 
B = E x v/c in the coordinate frame of the electron. This 
field partially polarizes the electron spin in the direction 
that is opposite for the electrons passing the nucleus on 
the right and on the left. Consequently, one can achieve 
a spatial separation of spin-up and spin-down electrons 
when unpolarized electron beam passes through unpolar- 
ized atomic target. 

A number of microscopic models have been developed 
that extended these ideas to solids. Various "extrinsic" 
(due to impurities) and "intrinsic" (impurity-free) mech- 
anisms have been studied with the aim to explain quan- 
titatively the spin Hall effect in non-magnetic conductors 
, as well as to explain the anomalous Hall effect in mag- 
netic materials [1] . Most of the existing theoretical mod- 
els are based upon the Boltzmann-type kinetic equation 
that describes spin and charge transport under certain 
assumptions about the collision integral. There is also a 



vast amount of numerical work of disordered systems in 
lattices of finite size. So far these works have not offered 
any universal description of the spin Hall effect. Instead a 
variety of different spin Hall effects (spin precession, side 
jump, skew scattering) that are specific to the model of 
spin-orbit interaction, nature of scatterers, band struc- 
ture, boundary effects, etc. has been proposed. A simple 
single-electron picture of the intrinsic spin Hall effect, 
similar to the treatment of the conventional Hall effect 
within Drude model, has been suggested in Ref. |9|. In this 
picture the Mott scattering of the transport current by 
the crystal field appears naturally within the framework 
of the Aharonov-Casher effect The parameter-free 

expression for the spin Hall conductivity was obtained 
for cubic crystals in good agreement with experiments. 
In this paper we will extend the treatment of Ref. to 
non-cubic crystals. Straightforward experiments will be 
suggested that can test our predictions. 

Our approach is based upon the general form of the 
one-electron Hamiltonian that contains spin-orbit inter- 
action to 1/c 2 fllj ]: 
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It is exactly this Hamiltonian, with U(r) being the elec- 
tric potential of the atom, that is responsible for the Mott 
scattering of electrons in the atomic physics. In a solid, 
U(r) is the electrostatic crystal potential felt by a charge 
carrier. For certainty we will speak about electrons but 
the model will equally apply to holes. With an accuracy 
to 1/c 2 Hamiltonian |lj is mathematically equivalent to 
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H — (p--A, T ) +U(r) 
where e < is the charge of electron and 
A, = -— — (cr x VU) . 
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(We have used the mathematical fact that the action of 
the operator p on A CT is zero, eijkO'j'S/iVkU — 0.) Con- 
sequently, the orbital motion of electrons is affected by 
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the fictitious spin-dependent magnetic field: 



V x A CT 
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[V x (er x VU)] . (4) 



This fictitious field produces the same effect on the or- 
bital motion of electrons as the real magnetic field does 
in the conventional Hall effect, but with the Hall cur- 
rents having opposite directions for electrons with op- 
posite spin polarizations. The spin-dependent Lorentz 
force, F^o- = — (e/c)(v x Bo-), gives rise to the fictitious 
spin-dependent Hall electric field, 



E_ff(T — Rh (Bo- x j ) 



(5) 



where jo = envo is the transport current expressed 
through the concentration n and drift velocity vo of the 
electrons, and Rh — —l/(nec) is the Hall constant. If 
the electrons are polarized, the spin average of Eq. (j4j 
gives rise to the effective magnetic field, 
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and to the measurable Hall electric field: 
Eff = Rh (B e // x j ) . 
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In these formulas < £ < 1 is the polarization of the 
electrons and (...) denotes the space average. If polarized 
electrons are injected from, e.g., a magnetic metal, jo 
describes the flow of the injected electrons. 

In the presence of the real magnetic field B, Eq. @ 
should be replaced by 
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U(r) + ^ B *-B, (8) 



where A=(Bx r )/2, g is the gyromagnetic factor, and 
\xb is the Bohr magneton. Since the definition of A CT , Eq. 
(J3j) , already contains 1 /c, the cross-term proportional to 
A • Ao in Eq. © has the order 1/c 3 . Within the non- 
relativistic approximation of Eq. JT]), that has accuracy to 
1/c 2 , such a cross-term must be omitted. Consequently, 
in the non-relativistic theory the conventional Hall effect 
and the spin Hall effect are totally independent. For 
that reason and in order to emphasize the consequences 
of the spin Hall effect, we are considering below the case 
of B = 0. 

When writing down Eq. ([5|) and Eq. we made an as- 
sumption that (v x Bo) = (v) x (B CT ). Some justification 
of this assumption follows from the fact that the trajec- 
tory of the charge carrier does not correlate strongly with 
the quadrupole component of the crystal electric field 
contained in the expression for B CT . Another argument 
is based upon symmetry. Indeed, the only reason for 
(F^o) to be different from zero would be v = (v) / 0. 
Consequently, (Fl<t) should be first order on vo- Being 
perpendicular to the velocity, the force F l g does not do 



mechanical work on the charge. Neither should (F^o) 
with respect to the drift motion of the charges, render- 
ing the form (Fj>) = — (e/c)(v x B). It is natural to 
identify B with B e/ / of Eq. ©. 

To show that space averaging in Eq. © produces 
a non-zero result, one needs to compute (VjVjJ7(r)). 
Experiments performed to date have been done in cu- 
bic semiconductors and in aluminum that is also cubic 
[E 0, 01 ■ For a cubic lattice 



(ViVjU) = Ch 
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due to the cubic symmetry alone, with C being a con- 
stant. This constant can be found from the Laplace equa- 
tion: 
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where p(r) is the charge density that creates U(r). To 
make the right choice for p(r), we notice that the spin- 
orbit interaction becomes larger as the electron passes 
closer to the nucleus. Similar to the Mott scattering by 
individual atoms, the electric neutrality of the crystal as 
a whole [1 31 ] is irrelevant for our problem; the distances 
that matter are the ones where the screening of the elec- 
tric charge of the nuclei is not complete [14| . Conse- 
quently, the spatial average in Eq. ([6]) must be over short 
distances. Due to the periodicity of the crystal it can 
be computed over the unit cell. In that sense our sym- 
metry argument for the cubic lattice and other lattices 
studied below is similar to the argument used to com- 
pute the crystal field (magnetocrystalline anisotropy) in 
magnetically ordered crystals [HI]. 

In line with the conventional approach to solids 
we choose U (r) as the potential formed by a cubic lattice 
of ions of charge — Ze > 0. Then (p) = —Zeno = —en 
where no and n = Zno are concentrations of ions and con- 
duction electrons, respectively. This gives C = 47re 2 n/3 
Q. Those who find this argument too simplistic may 
want to compare it with the approach developed by 
Hirsch [16]. In a model that replaces moving spins with 
stationary electric dipoles Hirsch computed the same av- 
erage over a cubic lattice of charges numerically. To nine 
decimal places his result coincides with ours up to a factor 
of 2 that can be traced to the difference in the expression 
for the fictitious magnetic field 17|. Note that in a mi- 
croscopic model the right-hand side of Eq. (fT0|) contains a 
sum over delta-functions. It is therefore likely that at the 
microscopic level the spin Hall effect originates from the 
singularity of the Coulomb potential. Since our model 
does not treat the screening effects rigorously and does 
not take into account the interaction between the elec- 
trons, our result for C can only be valid up to a factor of 
order unity. 

Substitution of Eq. §§§ into Eq. ^ gives 
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which provides the Hall field 
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Note that this formula does not contain any dependence 
on the concentration of charge carriers. Only the knowl- 
edge of £ and jo are needed to compute the Hall elec- 
tric field. This allows easy analysis of experiments with 
polarized electrons in cubic conductors. One can see a 
potential problem, however, with interpreting the Hall 
voltage in Eq. (fT2|) as a spin Hall effect. Indeed, the ef- 
fective field in Eq. (|11|) equals the magnetic field that 
the polarized electrons would produce in a spherically 
shaped body; with 47r/3 being the demagnetizing factor. 
The magnetic field due to polarization of electron spins 
in the actual sample should be different by a factor of 
order unity. Nevertheless, given the experimental uncer- 
tainties, it may be difficult to practically distinguish the 
spin Hall effect described by Eq. (fT2"]l from the conven- 
tional Hall effect due to the magnetic field produced by 
the electron polarization. 

In fact, the above-mentioned controversy is a conse- 
quence of the cubic symmetry. Indeed, in the conven- 
tional Hall effect only the orbital motion of the electron 
matters, the electron spin is irrelevant. On the contrary, 
the spin Hall effect arises from the spin-orbit term in the 
Hamiltonian. Consequently, there must be a clear way 
to distinguish between the two effects. As we shall see 
below, the non-cubic crystals present such a possibility. 
Consider, e.g., a tetragonal crystal with n being the unit 
vector in the direction of the c-axis. By symmetry, Eq. 
© should be now replaced with (a = b) 



(ViVjU) = C[k a §ij + (k c - k a )mnj 
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where k a ^ k c are factors of order unity. Working out 
the cross-products in Eq. ^ one obtains 
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'{ka + k c )£ + (k a - fc c )(n • £)n] . (14) 



The new feature is the component of the effective field 
along the tetragonal axis. (Notice the analogy with the 
magnetic anisotropy field in a magnetically ordered uni- 
axial crystal 151].) 

Substitution of Eq. (dU into Eq. gives 
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This formula must be also correct for a hexagonal crystal, 
with n along the hexagonal axis. Its remarkable property 
is that due to the second term the Hall voltage can be 
produced in the sample even when the electrons are po- 
larized along the direction of the transport current. This 
would be a clear manifestation of the Hall effect due to 
the spin, in contrast with the ordinary Hall effect due to 



the conventional spin-independent Lorentz force on the 
transport current. One possible geometry of the experi- 
ment is shown in Fig. [1] 




FIG. 1: Color online: Geometry of the proposed experiment. 
The crystal anisotropy axis is at an angle with the film. Due 
to the spin Hall effect the Hall voltage appears even when the 
electron spins are polarized along the transport current. 



Until now we have studied the spin Hall effect in a 
system of partially polarized electrons. Alternatively, 
one can study the spin currents that would polarize the 
boundaries of the sample in the absence of the polariza- 
tion in the bulk. In fact, this is how the spin Hall effect 
was initially observed P, Q • The spin current is defined 
[7j as the one-particle expectation value of 
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In general it is not conserved. However, due to the rel- 
ativistic smallness of spin-orbit interaction the reversal 
of the electron spin occurs only in a small fraction of 
scattering events. Consequently, in a small sample at 
low temperature the charge carriers typically reach the 
boundary of the sample before scattering reverses their 
spins. In this case Eq. (fTrj|) provides a useful concept for 
the study of the spin accumulation at the boundaries. 

Since Eq. (fl"6|) contains Pauli matrices, the non-zero 
expectation value of jik is provided by the spin-dependent 
part of the velocity operator. The latter is given by 



er 
m 



(17) 



where r is the scattering time. Consider, e.g., an or- 
thorhombic crystal. By symmetry the principal axes 
of the second-rank tensor (VjVjt/) should be directed 
along the a, b, c axes of the crystal. Choosing the axes of 
the coordinate frame along the crystal axes, we present 
(Vi'VjU) in a diagonal form: 



I] 1 



(18) 



with k x ^ y _ z being generally unequal factors of order unity. 
Substitution of Eq. (JTSJ) into Eq. Q then gives 
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Finally, with the help of equations (O, (fT6|) , (fT7)) . and 
(fl~9l). one obtains 
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where tro = e 2 nr jm is the usual charge conductivity. Ac- 
cording to this equation, the strength of the spin current 
as compared to the charge current is determined by the 
factor hao / (mc 2 ) which also determines the ratio of spin 
Hall and charge conductivities For good metals at 
low temperature this ratio can be of order 10~ 4 . 

For a cubic crystal k x = k y = k z = 1 and Eq. (j2"0|) 
reduces to the expression 
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which is independent of the orientation of the crystal 
lattice with respect to the transport current. This may 
present a problem for distinguishing the polarization of 
the sample boundaries generated by the spin Hall effect 
from the spin polarization generated by the Zeeman effect 
due to the magnetic field of the current. This controversy 
can be resolved by studying the spin Hall effect in a non- 
cubic crystal. Consider, e.g., two geometrically identical 
conducting strips, one cut along the aa-plane and the 
other cut along the ac-plane of a tetragonal crystal, with 
the z-axis being perpendicular to the plane of the strip 
and the transport current being along the y-axis. It is 
easy to see from Eq. ([20| that the spin current j zx = 
[nhao / (3mc 2 )](k x + k y )jo, describing the flow of a z along 
the x-axis, is different for the two strips. The ratio of 
these spin currents for the same value of the transport 
current is given by j zx {aa)/ j zx (ac) = 2k a /(k a + k c ) ^ 1. 
Consequently, the spin polarizations of the boundaries 
generated by the spin Hall effect will also be different 
for the two strips, while polarizations generated by the 
Zeeman effect due to the magnetic field of the transport 
current will be the same. 

In Conclusion, we have studied the intrinsic spin Hall 
effect in non-cubic crystals. The unique dependence of 
the effect on the crystal symmetry permits geometry of 
experiment in which the spin Hall effect can be unam- 
biguously distinguished from the effects caused by the 
orbital motion of charge carriers and by the magnetic 
field of the transport current. 
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